Abstract. We have derived the galaxy luminosity func- 
tion (GLF) in the cluster of galaxies Abell 496 from a 
wide field image in the I band. A single Schechter func- 
tion reproduces quite well the GLF in the 17< Iab <22 
(—19.5 < Mi < —14.5) magnitude interval, and the power 
law index of this function is found to be somewhat steeper 
in the outer regions than in the inner regions. This result 
agrees with the idea that faint galaxies are more abundant 
in the outer regions of clusters, while in the denser inner 
regions they have partly been accreted by larger galaxies 
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1. Introduction 

Galaxy luminosity functions (hereafter GLF) are funda- 
mental to analyse the properties of galaxies in clusters. In 
a number of cases, it is impossible to fit the entire GLF 
with a single Schechter function: there appear to be two 
components in the GLF, one for the bright galaxies - a 
gaussian distribution, and another for fainter galaxies - 
a power law or a Schechter function (see e.g. Godwin & 
Peach 1977, Biviano et al. 1995, Durret et al. 1999a). This 
suggests that there are at least two populations of galaxies 
in clusters, which do not vary strongly from one cluster to 
another, since the dip between both curves falls roughly 
at the same absolute magnitude in several clusters (see 
e.g. Tab. 2 in Durret et al. 1999a). Besides, at faint mag- 
nitudes, the slope of the GLF can be steeper in the out- 
skirts of clusters i.e. in less dense environments, and flat- 
ter near the cluster center (Lobo et al. 1997, Driver et al. 
1998, Adami et al. 1998, 2000, Andreon 2002, Beijersber- 
gen et al. 2002). This can be interpreted as due to the fact 
that in dense environments, small (and faint) galaxies are 
more likely to be accreted by larger ones. Moreover, they 
have also probably suffered repeated tidal interactions on 
their way towards the cluster center, consequently being 
dimmed or even disrupted (see e.g. Moore et al. 1996, 
Phillipps et al 1998, Kajisawa et al. 2000). 

We have performed a first analysis of the GLF of Abell 
496 (Durret et al. 2000) and intend to reobserve Abell 496 
spectroscopically with the VLT and VIRMOS; as a prepa- 
ration, we asked R. Ibata and C. Pichon to obtain for us 
a wide image of this cluster with the CFH12K camera at 
CFHT. We present below our analysis of the GLF in dif- 
ferent regions at various distances from the cluster center. 

2. The data 

2.1. Observations, reduction and detections 

Abell 496 (at redshift z=0.033, giving a distance modulus 
of 36.5, assuming H = 50 km s~ 4 Mpc~ 4 and q = 0, 



1 1 1 1 1 1 1 1 


1 




i 1 i i i 1 
ft ' * 






■ ■ 

































-4000 2000 2000 4000 

X (arcsec) 

Fig. 1. Positions of all the galaxies (relatively to the cD) 
with redshifts in Abell 496 from the Durret et al. (1999b) 
catalogue (crosses). The large rectangle shows the size of 
the present image; the small squares forming a sort of cross 
in the center correspond to the CCD catalogue by Slezak 
et al. (1999); the squares drawn with dotted lines indicate 
the fields covered by Molinari et al. (1998). 



as also used throughout this paper) was observed at the 
CFHT with the CFH12K camera in the Mould I band on 
February 20, 2001. Two images of 5 minutes exposure time 
each were obtained. The CFH12K camera is made of 12 
2Kx4K CCDs (hereafter, we will label the various CCDs 
from A to L anticlockwise from the south east corner). 
A global sketch of our field, together with those previ- 
ously observed by Molinari et al. (1998) and by our team 
are shown in Fig. [l], superimposed on the positions of the 
galaxies with redshifts belonging to the cluster (see the 
Durret et al. (1999b) catalogue). The pixel size of our im- 
age is 0.206 arcsec and the seeing 0.75 arcsec. The inter- 
ference fringes were corrected for and the photometrical 
calibration was estimated from the observation of the Se- 
lected Area 101 in the I Kron-Cousins system (Landolt 
1992), then converted to the Iab system by Lab =1+0.456 
(Fukugita et al. 1995). All our I magnitudes will here- 
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after be Iab magnitudes. The images were co-added and 
checked astrometrically at the TERAPIX data processing 
center, leading to a final image of 12365x8143 pixels, or 
42.45x27.96=1187 arcmin 2 in the East-West and North- 
South directions respectively. 

The sources were extracted using the SExtractor pack- 
age (Bcrtin & Arnouts 1996). Saturated objects (SEx- 
tractor flag > 4) were eliminated. The total number of 
objects thus eliminated was 154 in the magnitude range 
17< Iab _22 of interest here (see below). Their number 
varies from to 19 objects from one CCD to another, ex- 
cept for CCD L which has 69 (a first reason to discard this 
CCD). So, CCD L apart, these numbers are at most 4% 
of the total number of galaxies used to derive the GLF in 
each CCD and therefore eliminating them cannot strongly 
influence our results. 

In order to avoid false detections at the edges of each 
of the 12 CCDs, the catalogue of detected objects was 
limited to an area decreased by 15 pixels on all sides of 
each CCD. We thus obtained a final catalogue of 37058 
objects in a total area of 1123 arcmin 2 (or 33596 objects 
in 1031 arcmin 2 if CCD L is excluded). This catalogue 
will be made available in electronic form, with the follow- 
ing Cols.: (1) running number; (2) and (3) right ascen- 
sion and declination; (4) major axis a in arcsec; (5) major 
axis position angle; (6) ellipticity; (7)-(8) integrated ellip- 
tical Kron Iab magnitude and corresponding error; (9)- 
(12) aperture magnitudes within 15, 10, 5 and 3.64 pixels 
respectively (3.64 pix=0.75 arcsec, the FWHM of the see- 
ing). All these parameters are those estimated with SEx- 
tractor. Since we have only one filter and a rather short 
exposure time we cannot give any accurate morphologi- 
cal information. However, the combination of some of the 
information provided in this catalogue, such as the total 
magnitude versus the magnitude in an aperture having a 
diameter equal to the seeing FWHM, can be used to ob- 
tain a first order morphological indication in the form of 
a concentration parameter. 

We have checked our photometry with data in the lit- 
erature. First, we identified 19 bright galaxies in common 
with the LEDA data base ( http: / /lcda.univ-lyonl ,fi ) . We 
find a mean value < Bt — Iab >= 1-89 (cr=0.36). As- 
suming Iab=I+0.456 this gives < Bt — I >= 2.35, in 
agreement with the value given for elliptical galaxies by 
Fukugita et al. (1995): B-I=2.27. Second, we retrieved the 
Moretti et al. (1999) catalogue for Abell 496 (which is 
broader and deeper than our previous R band catalogue) 
in the Simbad data base. For 36 galaxies we obtain <r- 
Iab >=0.70 (dispersion 0.53) and <i-Iab >= 0.55 (dis- 
persion 0.55); these values correspond to <r-I>=1.16 and 
< i — I >= 1.01 with the above conversion, to be com- 
pared with the respective values of 1.04 and 0.75 given by 
Fukugita et al. (1995). Therefore, the agreement with the 
magnitudes of the Moretti catalogue is correct, despite a 
possible shift by at most ~ 0.25 magnitudes. The agree- 
ment is good with the LEDA catalogue. 



Fig. 2. Variation of the completeness level as a function 
of I magnitude for 24 different regions of the image (each 
CCD was split horizontally into 2 equal sub-areas labelled 
t (top) and b (bottom)). The thick curve shows the com- 
pleteness averaged over all CCDs but L. The thick hori- 
zontal line indicates the 90% completeness level. 



2.2. Completeness 

We have carried out simulations to compute the complete- 
ness level. We added artificial objects (similar to real ob- 
jects) to our image and measured the fraction of these 
objects recovered by the SExtractor package as a function 
of magnitude and location in the image. For this, we used 
a code created by J.M. Deltorn and already applied to 
the CFDF survey (e.g. McCracken et al. 2001) to com- 
pute the star detection completeness level. We modified 
this code in order to have a more realistic representation 
for galaxies, and used a gaussian profile with a FWHM of 
3 times the mean seeing of our observation. The results 
are given in Fig. |. This fi gure represents the percentage 
of completeness level as a function of I magnitude for 24 
different regions of the total image (each CCD was split 
horizontally into 2 equal sub-areas). CCD L is significantly 
less complete than the other CCDs, because of many dead 
columns, and was removed in the following analyses. For 
the other CCDs, the mean completeness level is close to 
90% up to Iab =22, which will be taken as the 90% com- 
pleteness limit for our catalogue. We also obtained another 
estimation of our completeness level by comparing our ob- 
servations with those of the CFDF survey (e.g. McCracken 
et al. 2001), which used the same instrument and filter, 
but with an exposure time of 5.5 hours and a seeing of 1 
arcsec, and was complete up to 1=25.6. Their completeness 
limit rescaled to our exposure time gives a completeness 
level at 1=21.8, in agreement with our simulations. We 
will not attempt to correct our counts for incompleteness 
at magnitudes fainter than Iab =22, and we will hereafter 
limit our analysis to Iab _22. 

In order to estimate the influence of "crowding" we 
computed the number of galaxies susceptible to be masked 
by bright galaxies. For 18< Iab _22 there are 9726 galax- 
ies in an area of 1.007 x 10 8 pixels 2 (entire field). As a 
conservative approach, we consider that galaxies brighter 
than Iab = 18 and with surfaces larger than 600 pixels 2 



4 



Durret et al.: The galaxy luminosity function in the cluster of galaxies Abell 496 



ckground galaxy counts 



Postman galaxy counts 

A496 outer zone 



Fig. 3. Field galaxy counts from the VIRMOS survey 
(dashed line) and from Postman et al. (1998) (full line) a 
magnitude shift of +0.456, normalized to the total area of 
1123 arcmin 2 covered by our catalogue. The galaxy counts 
derived from the outer regions of our Abell 496 field as ex- 
plained in the text (Sect. 3.2.) are also shown. 



can mask faint galaxies. The total surface covered by these 
galaxies is 24130 pixels 2 , leading to a number of galaxies in 
the 18< Iab <22 magnitude interval that can be masked 
of the order of a few. Therefore, the influence of crowding 
on our study appears to be negligible. 

3. Estimating the background contamination 

Since the galaxy-star separation becomes difficult for mag- 
nitudes Iab > 20, we decided to subtract the star and 
background galaxy contaminations statistically. 

3.1. Star counts 

In order to subtract the stellar contribution from our 
Galaxy, we produced a catalogue of stars using the Be- 
sancon model (Gazelle et al. 1995) in the I Kron-Cousins 
band in the direction of Abell 496. The relation between 
the magnitudes measured with the two filters is again: 
lAB=lBcsancon + 0.456. The uncertainty on the star counts 
is smaller than 10% for Iab <22 (A. Robin, private com- 
munication), and the contribution of stars remains small 
in any case (less than one tenth of the total counts in 
our image, all objects considered), so star counts cannot 
influence our galaxy counts by more than a few percent 
in the magnitude interval 17< Iab < 22. An alternative 
way to perform this correction would have been to use the 
VIRMOS star counts, which have the advantage of having 
been obtained with the same CFH12K camera and filter, 
but they are not representative of the stellar counts in the 
direction of Abell 496, due to their lower Galactic latitude. 



as far as possible from the cluster center, subtract to it 
the star contribution and subtract the resulting galaxy 
counts to our data. We have extracted such an "outer 
zone" by putting together the data of the left half of CCD 
A and the right halves of CCDs F and G, and derived the 
galaxy counts in this zone. As can be seen in Fig. |^, these 
counts are higher than those issued from two independent 
field surveys (see below), once all are normalized to the 
same surface area (the total size of our image, i.e. 1123 
arcmin 2 ), suggesting that the Abell 496 "outer counts" 
thus produced still contain a significant fraction of cluster 
member galaxies, as confirmed by the positions of cluster 
galaxies in Fig. |l[ So we would obviously overestimate the 
background if we took it in this "outer zone" . Note that 
our image covers a total region of 2.404x 1.584 Mpc 2 , while 
the r2oo radius calculated as in Carlberg et al. (1997) with 
a velocity dispersion of 715 km/s (Durret et al. 2000) is 
2.5 Mpc, confirming that the cluster contribution in the 
outer regions of our image is non negligible. Besides, the 
fact that the last point of the counts in the outer zone 
(Iab =23.75) does not merge with any of the field survey 
background counts described below adds still another rea- 
son to reject this type of background subtraction. 

We therefore decided to subtract the background 
contribution taken from the VIRMOS survey galaxy 
counts (McCracken et al. in preparation). However, the 
background galaxy contribution taken from this survey 
shows an excess of objects at magnitudes brighter than 
Iab=18.5, and cannot be directly subtracted to our counts 
either (see Fig. ||). The comparison of the galaxy counts 
by Postman et al. (1998) to the VIRMOS galaxy counts 
shows similar slopes for Iab > 18.5 (see Fig. ||), with 
a magnitude shift due to the fact that the VIRMOS 
counts are in Iab magnitudes while the Postman counts 
(which match well other counts such as those of Ca- 
banac et al. 2000) are in Cousins I magnitudes. Shifting 
the Postman counts by 0.456 magnitude (see Sect. 2.1) 
gives a good agreement between both background counts 
for Iab >18.5. We will therefore subtract to our data 
the Postman galaxy counts shifted by 0.456 magnitude 
for 17< Iab <18.5 and the VIRMOS galaxy counts for 
18. 5< Iab <22. The star and background galaxy subtrac- 
tions were performed in bins of 0.5 magnitude. The error 
bars that we indicate in the plots are simply computed 
as the square root of the number of galaxies in the corre- 
sponding magnitude bin. 



4. Fitting method and results 

4-.1. Fitting method 



3.2. Field galaxy counts 

An obvious way to correct for the background galaxy con- 
tribution would be to extract from our image a region 



Schechter function fits were performed using an IDL code 
based on the curvefit function, which uses a gradient- 
expansion algorithm to compute a non-linear least squares 
fit to a given function; this routine gives the best fit pa- 
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Galaxies towards A496 
Cabanac + Virmos 
Galaxies in A496 




22 24 



Fig. 4. Galaxy counts in the direction of Abell 496 (dot- 
dashed line), background galaxy counts - see Sect. 3.2 
(dashed line), and final galaxy luminosity function for 
Abell 496 after background subtraction (full line) with its 
error bars. 

rameters and respective errors of the Schechter function: 

$(M)dM = a-$*10°- 4 ( m *- m » q+1 ' exp(-10 a4 ( M *- M >) 

where <&* is the normalisation, M* the characteristic ap- 
parent magnitude, a the slope of the faint end of the lumi- 
nosity function and M the apparent magnitude of a given 
galaxy in the I band. 

4-2. Overall galaxy luminosity junction 

The resulting GLF for the entire field (after eliminating 
CCD L, the top half of CCD B: Bt and the bottom half of 
CCD D: Db) is shown in Fig. |j. A single Schechter func- 
tion was fit in the same magnitude interval that we have 
been using (17< Iab <22, -19.5 < Mi < -14.5), giving a 
power law index slope a = —1.79 ±0.01. Results are given 
111 Table | for all regions that we explore. For each case, 
we indicate the reduced x 2 ; all values are about 1 or lower, 
indicating that the fits are correct. Since M* is brighter 
than the lower limit of the magnitude interval considered 
here, it is not well constrained and we will therefore focus 
our discussion on the values of the slope a only. The cor- 
relation between a and M* is shown as confidence ellipses 
in Fig. H for various regions. The ellipses confirm that the 
error bars that we give on a in Table | are realistic. Note 
that the error bars in Table 1 are la error bars and corre- 
spond to the innermost ellipses in Fig. |[ We made tests 
on region CDIJ (see below), to see how an underestimate 
of the error bars could modify the slope of the GLF. For 
this, we multiplied the error bars by factors of 2 and 10 
and found that a remains unchanged (even though our 
fitting procedure takes the error bars into account), while 
the uncertainty on a increases to ±0.05 and ±0.1 respec- 
tively, instead of the previous value of ±0.03. In order for 
our results to lose significance, we would have to multi- 
ply the error bars by 10, which seems an unrealistically 
large number (as seen for example from the scatter in the 



Region 


Nb. 
gal. 


Area 
(arcmin 2 ) 


a 


/Cred 


All* 


4052 


1031 


-1.79 ±0.01 


7.26/7 


ABK 


1217 


281 


-1.93 ±0.02 


4.14/7 


CDIJ 


1542 


378 


-1.75 ±0.03 


1.50/7 


EFGH 


1652 


372 


-1.98 ±0.03 


1.80/7 


CenL 


426 


94 


-1.82 ±0.04 


1.80/7 


CenM 


292 


53 


-1.79 ±0.07 


0.48/7 


CenS 


129 


24 


-1.60 ±0.25 


1.50/7 


A 


409 


93 


-2.05 ±0.05 


3.60/7 


B 


402 


94 


-1.73 ±0.04 


3.18/7 


C 


312 


95 


-1.73 ±0.05 


3.60/7 


D 


353 


94 


-1.87 ±0.04 


3.36/7 


E 


364 


94 


-1.68 ±0.13 


1.98/7 


F 


460 


92 


-2.03 ±0.05 


2.34/7 


G 


297 


92 


-1.87±0.11 


2.76/7 


H 


534 


94 


-2.03 ±0.05 


4.08/7 


I 


374 


94 


-1.60 ±0.10 


0.66/7 


J 


505 


94 


-1.78 ±0.08 


2.46/7 


K 


407 


94 


-1.94 ±0.05 


4.08/7 



All but L,Bt,Db 



196 GT,F 



'_' 1 1 1 1 1 1 1 1 1 1 1 '_ 


<_< \ < < < \ < < < \ < < <_ 

R 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

15 

E 


ARK lE^jSf 

- 7^1 

1 , , , 1 , , , 1 , , , 1 


G 

_ H 


J 

" K _---\ 



IB 20 22 24 1 a T20 22 24 18 20 22 24 
IAB 

Fig. 6. Galaxy luminosity function in various regions of 
Abell 496 (see text and Table 1). 



Metcalfe et al. Fig. 13). A factor of 2 seems much more 
probable, and in this case the difference in slopes between 
CDIJ and other regions remains significant. 

4-3. Galaxy luminosity function in three large regions 

We then divided the cluster into three regions of roughly 
equal surface, CDIJ surrounding the cluster center (CCDs 
C, D, I and J), ABK (CCDs A, B, K) towards the East and 
EFGH (CCDs E, F, G, H) towards the West. The GLFs 
in these three regions are shown in Fig. ||. The fit of the 
GLF by a Schechter function in region CDIJ is displayed 
in Fig. [7] (the fit was done in the interval 17< Iab <22 
even though the figure shows the GLF in a larger range of 
magnitudes, for which we also extrapolated this best fit). 

The slope of the GLF is found to be flatter in region 
CDIJ, the Schechter law slopes being a = —1.75 ± 0.03, 
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Fig. 5. Correlation between the two parameters M* and a of the Schechter function for six regions: all the cluster, 
ABK, CDIJ (top row), EFGH, F and C (bottom row). The isocontours are 1, 2 and 3a respectively. 
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Fig. 7. GLF in the central region (CDIJ) of Abell 496 with 
the best fit Schechter function superimposed. 

-1.93 ± 0.02, and -1.98 ± 0.03, for regions CDIJ, ABK 
and EFGH respectively (see Table Q). 

4- 4- Mapping the parameters of the galaxy luminosity 
function 

The GLFs in the 11 CCDs (L excluded) are shown in Fig. || 
and Table |l|. Here also, the Schechter function slope may 
be steeper in the outer regions of the cluster, but better 
statistics are obviously required. 

Finally, we selected three rectangular concentric re- 
gions of different sizes centered on the intersection of 
CCDs C, D, I and J (or, roughly the position of the cD 
cluster galaxy). CenL (Large), is of the size of a CCD 
and covers one quarter of CCDs C, D, I and J. CenM 
(Medium), covers 9/16 of the area of CenL. CenS (Small), 
covers 1/4 of CenL. The Schechter fits of the luminos- 
ity functions give slopes a = -1.82 ± 0.04, -1.79 ± 0.07 
and —1.60 ±0.25 for CenL, CenM and CenS respectively, 
consistent with that in the CDIJ area within error bars. 
Region CenS may exhibit a flatter slope, but this needs 
confirmation since the uncertainty is very large. 

5. Discussion and conclusions 

We have derived the GLF in various regions of Abell 496. 
The slope of the Schechter function fit is always found to 
be steep (between —1.60 and —2.05). Since such a steep 



slope could be due to several artefacts, we will discuss the 
validity of our results. First, the background counts could 
have been underestimated. However, the good agreement 
of the various background counts (VIRMOS, Postman, 
Cabanac), and the fact that the subtraction is mainly that 
of the VIRMOS counts (in the interval 18. 5< I A b <22), 
made with the same instrument, filter and magnitude sys- 
tem as ours, tends to suggest that this is not the case. 
Second, the number of faint galaxies may have been over- 
estimated; for example, we may have confused globular 
clusters with galaxies at faint magnitudes, as explained in 
detail by Andreon & Cuillandre (2002). However, we limit 
our sample to Iab=22, where such effects should not be 
too strong. Third, our Iab magnitudes may be too bright 
by 0.25 magnitude, as suggested by the difference with the 
Moretti et al. data (see Sect. 2.1). We tried to fit the GLF 
in several regions after shifting the Iab magnitudes by 
0.25 (before subtracting the background) and find slopes 
a = -1.68 ± 0.05, -1.85 ± 0.03 and -1.94 ± 0.04 for re- 
gions CDIJ, ABK and EFGH respectively, instead of the 
previous values of —1.75,-1.93 and —1.98. Therefore, al- 
though the values change a little, the slope remains flatter 
in CDIJ. 

We then used our Iab catalogue limited to the region in 
common with Molinari et al. (1998) and made a Schechter 
fit as described above. The Molinari et al. zone partially 
covers our CCDs J,I,C,D,E and F, with a main concen- 
tration towards CCD E. A Schechter fit in the same mag- 
nitude interval gives a slope a = —1.71 ± 0.06, close to 
the value of -1.68 ± 0.13 found in CCD E. A shift of I AB 
by 0.25 magnitude as above gives a = —1.68 ± 0.09, in 
perfect agreement with Molinari. Fits in broader magni- 
tude intervals give: a = —1.66 ± 0.05, —1.64 ± 0.04 and 
— 1.57 ±0.03 for the magnitude intervals 17< Iab <22.5, 
17< Iab <23, and 17< Iab <23.5 respectively. Molinari 
et al. give a slope a — —1.49 ±0.04 in the i band, but men- 
tion that a magnitude correction allows them to reach a 
slope as steep asa = —2.0. We therefore believe that our 
results are consistent with theirs. Note that such a slope 
is not much steeper than found e.g. in Coma (Lobo et al. 
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1997) or in Abell 665 (De Propris et al. 1995). This could 
indicate an excess of faint red galaxies, but to ascertain 
this hypothesis it would be necessary to derive the GLF 
in Abell 496 in other filters from samples of comparable 
quality (covered area and depth). As still another test on 
the robustness of our results, we reanalyzed the GLF in 
the CDIJ region. For this, we reduced the number counts 
by 10, 20, 30 40 and 50% for I AB >20 (below I A b=20 the 
counts remained unchanged) and made fits of these new 
GLFs. Results are given in Table 2. They show that the 
slope changes strongly only if the counts are reduced by 
at least 30%, an unrealistic number. Besides, in order to 
account for the difference in slope of 0.2 that we observe 
for example between regions CDIJ and EFGH, we would 
need to make an unrealistically large error of 50% on the 
counts. We are therefore confident that both the absolute 
values of a and their variations from one zone to another 
are robust. 

Our second result is that the slope of the Schechter 
function fit tends to be steeper in the outer regions of the 
cluster, as already observed in other clusters (see refer- 
ences in Sect. [I]). Such a variation of a can be interpreted 
as due to the fact that faint galaxies are accreted by larger 
ones preferentially in the inner parts of clusters, where 
the galaxy density is higher, therefore inducing a lack of 
faint galaxies and a flattening of the GLF in the inner 
regions. Moreover, galaxies are likely to have suffered re- 
peated tidal interactions on their way towards the cluster 
center, consequently being dimmed or even disrupted in a 
scenario of galaxy harassment (Moore et al. 1996). 

The next step is obviously to confirm these results 
through deep multiband imaging and/or spectroscopy 
that would make the background subtraction more secure 
and would allow us to compare the GLFs in various filters. 
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